Reynolds Pressure and Relaxation in a Sheared Granular System 
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We describe experiments that probe the evolution of shear jammed states, occurring for packing 
fractions <j)s < <j> < 0j, for frictional granular disks, where above <j)j there are no stress- free static 
states. We use a novel shear apparatus that avoids the formation of inhomogeneities known as shear 
bands. This fixed cj) system exhibits coupling between the shear strain, 7, and the pressure, P, which 
we characterize by the 'Reynolds pressure', and a 'Reynolds coefficient', R(<j)) = (d 2 P/d'j 2 )/2. R 
depends only on 0, and diverges as R ~ (</> c — 4>) OL : where <j) c ~ 0j, and a c± —3.3. Under cyclic shear, 
this system evolves logarithmically slowly towards limit cycle dynamics, which we characterize in 
terms of pressure relaxation at cycle n: AP ~ — f3 ln(n/no). /3 depends only on the shear cycle 
amplitude, suggesting an activated process where /3 plays a temperature-like role. 
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Much recent work has focused on jamming and glassy 
behavior in disordered systems of particles, including 
granular materials, colloids, foams and molecular glass 
formers. Jamming, as encapsulated in the Liu-Nagel jam- 
ming diagram [IH3], has provided key perspectives. How- 
ever, Bi et al. [4j [5] recently showed that frictional ma- 
terials, e.g. most common granular systems, differ from 
the L-N scenario: shear strain, 7, can 'shear jam' [5] low 
density states. There is a lowest packing fraction <j)j such 
that below (above) this density, there are (no) zero-stress 
states [5]. Application of shear to a zero-stress state in 
0s < <j> < <t>j leads to highly anisotropic contact and 
force networks, and to non-zero shear stress, r, and pres- 
sure, P. Here, r = (<Ji — <T2)/2 and P — (a\ + cr 2 )/2, 
where the Gi are the principal stresses of the 2D stress 
tensor, a. Starting from zero stress, the system tra- 
verses a fragile regime, and with additional shear strain, 
the system arrives at a fully jammed state where the 
force/contact networks percolate in all directions. There 
are accompanying evolutions in the average contact num- 
ber per particle, Z, and the fraction of non-rattler parti- 
cles, fNRi where the latter provides a particularly clear 
signature of shear jamming. Improved understanding of 
shear jammed states is crucial for both a better under- 
standing of jamming for frictional materials, and to shed 
light on the complex rheology of dense granular flows [6] . 

A complication for sheared frictional materials is the 
formation of dilated localized shear bands, where most of 
the shear strain is confined. System-wide measures may 
tend to reflect the band properties rather than the whole 
system, making it difficult to interpret experiments. In 
this Letter, we show how to experimentally perform shear 
in a 2-D frictional granular material so as to avoid shear 
bands. We then use this method to characterize the me- 
chanical response and dynamics of shear jammed pack- 
ings. 

At the heart of shear jamming are classic studies by 



Reynolds, who showed that under fixed pressure, gran- 
ular systems can dilate in response to shear 0. For a 
system sheared at fixed density, as here, a related ef- 
fect occurs: the stresses a respond to the shear strain. 
We find that P increases roughly as 7 2 , which we de- 
scribe by a "Reynolds coefficient", R = (d 2 P / d^) /2. 
We find that R depends only on 0, and it provides a 
simple parametrization of the coupling between P and 
7. R seems to diverge as <j) approaches (j) c ~ </>j, thus 
identifying a special role for <\>j for the shear jamming 
states. 

An additional key observation from this work is that 
for (\>s < (\) < the stress response to cyclic shear strain 
shows slow relaxational dynamics to a limit cycle, that 
depends on driving. The deviation from a limit cycle, 
measured by pressure, shows a logarithmic decay over 
time/cycle number. The data for stress relaxation exhibit 
a simple scaling form, as developed below. 

Experimental Setup Key to these experiments is a novel 
apparatus that provides (simple) shear throughout the 
system, in contrast to wall-driven shear. The base of the 
apparatus consists of narrow, parallel, horizontal, and 
transparent slats. Shear is applied by deforming the slats 
and boundary uniformly in the 'y' direction, keeping the 
'x' dimension fixed at L, to provide uniform simple shear 
strain 7 = Ay/ L at constant packing fraction cj) (Fig.JlJi). 
On the slats rest ~ 1000 bi-disperse photo-elastic parti- 
cles (Vishay PSM-4) of diameters 12.7mm and 15.9mm; 
the slat width is of the order of the particle size. The rel- 
ative numbers of large to small particles is set to 1 : 3.3, 
in order to prevent crystallization. This bottom-assisted 
shear induces a linear shear profile, suppressing shear 
bands and the usual inhomogeneities. It is reminiscent 
of the SLLOD and related algorithms [8J for enforcing 
uniform shear in MD simulations. It bears some resem- 
blance to 3D experiments by Mueggenburg [9 , but with 
a key difference: in the Mueggenburg experiments, a slat 



2 




10 20 30 40 

x position/d 



FIG. 1: (Color online) (a) Setup schematics (b) The three 
close-up images that the camera captures at each step: parti- 
cle positions (upper), force response under polariscope (mid- 
dle), and particle orientation images under UV light (lower), 
(c) The x— and y- displacement of particles vs. their horizon- 
tal positions in the system, and (d) the coarse grained |15ill6| 
density profile after 27% linear shear. 




FIG. 2: (Color online) (a) Reynolds pressure P(7 2 ) observed 
in forward shear (see text) tests for cj) = 0.691 — 0.816. The 
data has a positive noise floor of ~ 0.5 N/m. (b) Reynolds 
coefficient R extracted from linear fitting, obtained from up 
to 54% forward shear (red squares), up to 27% forward shear 
(blue dots), and cyclic shear tests under limit cycle behavior 
(black triangles). The inset shows the same data on double 
logarithmic scales with (j) c = 0.841 ± 0.004. The error bar 
is smaller than the size of the symbols unless marked. The 
dashed line shows a fit to a power law. A line corresponding 
to an exponent -3.3 is also shown for reference. 



geometry was used, but the slat motion was not coordi- 
nated, and sustained uniform shear did not occur. 

The experiment is illuminated from below by circu- 
larly polarized uniform white light, and from above by a 
less intense UV light. A 22 Megapixel camera above the 
experiment records views with and without a circular po- 
larizer. We apply quasi-static shear strain in small steps. 
After each step, we pause and record three views of the 
system which respectively yield particle positions, photo- 
elastic responses, and rotations. Without crossed polar- 
izers, the edges of the particles are visible (Fig.[T]3 upper), 
and we use a circular Hough-transform technique [10] to 
determine particle centers with an accuracy of ~ 0.02<i. 
With a circular polarizer in front of the camera, we image 
the photo-elastic pattern of colored fringes within each 
particle, which encode the contact forces acting on each 
particle (Fig. [TJd middle). To determine the particle ori- 
entations, each disk is marked diametrically with a line 
of fluorescent dye, visible under UV light with the white 
light turned off (Fig. [IJ) lower). Changes in the bar ori- 
entations give particle rotations. The complete process 
of multiple strain steps, followed by imaging after each 
step, is fully automated, and we record up to 1000 steps 
of forward or cyclic shear per run. We extract the local 
particle stress by either a pattern-fitting approach [4] [11], 
yielding the complete contact network, particle forces, 
and stress tensor (e.g. P and r), or via G 2 , the local 
squared intensity gradient of the photo-elastic response, 
averaged on each particle [12j H3] , G 2 is a one-to-one 
function of P on the particle level, providing an efficient 
measure for P. For small/large data sets, we use the 
former/latter approach. 

Reynolds Effect As noted, a striking aspect of applying 
shear strain to a stress-free state for c^s < </> < </>j, is the 
generation of nonzero P and r, as in the shear jamming 
experiments of Bi et al. [5 . Here, we probe the evolution 



of shear jammed states by first forward-shearing the sys- 
tem. We prepared packings in a stress free initial state, 
for 0.691 < <j> < 0.816, where <j>j = 0.835 ± 0.005, and 
<ps ~ 0.75. We then quasi-statically shear the system by 
200 small strain steps of 0.27%, up to a total strain of 
7 = 54% [14] . We see most of the same general features 
of shear jamming, as reported by Bi et al.[5 , but here, 
particle tracking data (Fig. [l^) show that the shear is 
effectively linear and homogeneous across the entire sys- 
tem. Particle displacements and rotations relative to the 
uniform shear background are small. The locally coarse- 
grained density field [15j [16] (Fig. [l]i), shows no sign of 
a shear band or permanent inhomogeneities. 

For the larger cVs considered here, we could not ap- 
ply the full 54% strain because P became so large that 
the layer was unstable to out-of-plane buckling. If buck- 
ling occurred, we terminated the forward shear experi- 
ment. The forward shear results, Fig. [2^i, indicate that 
the shear-induced 'Reynolds pressure' increases roughly 
as 7 2 with a density dependent prefactor which we char- 
acterize by the 'Reynolds coefficient', 

R= (9 2 P/^)/2. (1) 

For linear isotropic elastic materials, no coupling between 
shear strain and pressure is expected. But, as we apply 
shear, the system becomes increasingly anisotropic, so a 
P— 7 coupling might be possible, as expressed by, dP/d^y. 
In our system, this derivative grows roughly as 7, and 
linear elasticity is not a particularly useful concept. R 
grows strongly with 0, apparently diverging at <p = <p c ~ 
0j, e.g. Fig. [2)3 and inset, which shows a log- log plot of R 
vs. A(j) = (j) c — (j>. A least squares fit to R = AA(j) a , yields 
a = -3.3 ±0.1 and (j) c = 0.841 ± 0.004. By contrast, <f> c 
lies in the range 0.83 < <j)j < 0.84, so here, <j) c is not 
distinguishable from 0j, which is also comparable to <pj 
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FIG. 3: (Color online) (a) P vs. 7 for a symmetric cyclic 
shear run with <j) — 0.825, which started from 7 = 0, and 
sheared between jmax = 2.25% and 2.25%. Only 

cycles 1, 2, 28 and 29 are shown in the plot, (b) r vs. 7 for the 
same run and the same shear cycles. (c)P vs. 7 at cycles 1, 
2, 28, 29 for a non-symmetric cyclic shear run (jmax = 4.5%, 
Jmin = 0) with the same density. 



for systems of frictionless 2D particles. For <p < 0.75, the 
system is very loose, and it does not form a percolating 
contact network, even after 54% strain. In this case, 
the very small R that we observe is likely due to the 
friction between particles and the base slats, and possibly, 
the finite size of the system. R((f)) behavior in this case 
deviates from the power-law behavior (Fig. [^(inset)). 
We identify (j>s — 0.75, the lower limit in this system for 
shear jamming. 

Limit Cycles To characterize the evolu- 
tion/reproducibility /relaxation of the stresses, we 
carried out multiple shear cycles. This also allowed us 
to determine R for <fi closer to 0j, where shear strains 
are limited due to buckling; we obtain good statistics by 
many smaller-amplitude strain cycles. The oscillatory 
shear experiments were started from initially stress-free 
states for 0's in the shear jamming regime, <ps < <t> < 4>J- 
In a cycle, we sheared by strain steps of 0.45% up to 
jmax in the 'forward direction', followed by a shear 
strain decrease (—0.45% per step) to a smaller strain, 
7mm- For symmetric shear cycles: j min = -jmax, and 
asymmetric shear cycles: ~f min ^ -Jmax- 

For symmetric cycles, P was symmetric about 7 = 0, 
approximately quadratic in 7, and virtually reproducible 
over many cycles, as shown in Fig. [3^i. The Reynolds co- 
efficient R(4>) followed the same trend as in the forward 
shear tests (Fig. J2Jd), further confirming the Reynolds 
effect. After transients, the shear stress r also fol- 
lowed a reproducible path over cycles, but unlike P, r 
was strongly hysteretic, with non-zero values at 7 = 
(Fig. [3]). There were 7's for which r = but P ^ 0, for 
example, in Fig. [3^1 and b, at 7 « 1%. However, in such 
cases, r coarse grained at smaller scales than the system 
size was locally non-zero, even though the global r was 



(e.g. because of spatial variations of the principal stress 
orientations). Due to length limitations, we consider only 
the dynamics exhibited by P, and we will present the full 
stress dynamics elsewhere. 

The evolution of P(j) for asymmetric shear cycles dif- 
fered from the symmetric case. Here, P(j) was initially 
asymmetric, but evolved towards a symmetric shape cen- 
tered around the mean strain, 7, after many cycles. Thus, 
the long term P — 7 dynamics was a limit cycle. The sys- 
tem relaxed quickly (slowly) to the limit cycle if sheared 
symmetrically (asymmetrically). Fig. [3^ shows an exam- 
ple of slow evolution, where a limit cycle was reached 
after about 28 cycles. In this case P("y) evolved to a 
symmetric shape, similar to the forward shear experi- 
ment, except for a vertical shift; i.e., the system did not 
reach a completely stress-free state at the mid-point of 
strain. However, a long-term limit cycle was still reached 
with the same Reynolds coefficient for the given strain 

amplitude, J A Jmax ~ Imin- 

Slow Relaxation For asymmetric strain cycles, 
AP(n) = P(jmax)-P(7min) was initially nonzero, but it 
decreased and ultimately vanished, within fluctuations, 
for n = no- When the limit cycle was reached, P was 
symmetric about 7 = (j m ax + 7mm) /2. The slow relax- 
ation of AP for asymmetric shear shows striking scaling 
behavior, which we characterize in terms of </>, 7 and 
the shear amplitude ja- Experiments to characterize 
this relaxation spanned 0's from above <ps to just below 
isotropic jamming <fij: 0.780 < <fi < 0.828, strain ampli- 
tudes of ja = 6.75,4.5,3,1.5% and a range of starting 
strains < 7 < 21.35%. Experiments were 100-500 cy- 
cles long; for convenience we measured G 2 only at 7 ma£C , 
jmin, and then converted G 2 to AP using a calibration. 
Fig. [4^i shows AP for a particular 7^4. 

For (j) in the shear jamming region, AP(n) decayed 
logarithmically slowly towards 0: 

AP(n)~-/31og(n/n ), (2) 

implying a natural 'time scale' for relaxation, no, that 
we obtained through least squares fits of the logarithmic 
part of the relaxation. All the relaxation data, for a 
given 7^, collapse onto a single curve when expressed in 
terms of n/no (Fig. [4)3), regardless of <j) and 7. The decay 
factor PipfAj differs for each 7^ (Fig.Efc), but AP//3 is a 
universal function of n/no, as in Fig.pp, which shows all 
~170 datasets. We emphasize the remarkable role that 
P(ja) plays, and the fact that it is independent of <p. 

We then consider what determines no- Eq. [2] implies: 
no = n • exp(AP(n)//3(7^)). Initially, at n = 1, AP = 
APo. According to the approximately quadratic relation 
between P and 7, APo is given by: APo = R{<t>){lmax ~ 
7min)/2 = R{4>)llA- Therefore, 

n = exp(R(</>)i ^ ). (3) 
Eq. [2] also implies an evolution dAP/dn = 
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stress ensemble, given the activated process character of 
the slow relaxation. Such a theory would need to explain 
some of the striking scaling properties observed here. In 
addition, we have not considered the properties of the 
shear stress under cyclic shearing, nor have we consid- 
ered the particle dynamics or details of the force/contact 
networks. We will present these results elsewhere. 
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FIG. 4: (Color online) (a) AP vs. number of shear cycles, 
for various runs with 7a = 2.7%, and 7 = 1.35% (solid line), 
5.85%(open circles), and 10.35%(open triangles). Different 
color indicates different <j). (b) AP vs. log (n/ no) for 3 different 
7a 's: 1.35%(open circles), 2.70%(open triangles), and 6.75% 
(solid lines). Data comes from runs with various <j> (shown 
by color) and various 7 (not shown), (c) The decay factor, 
/3, shows a strong increase with 7a- (d) The universal decay 
curve, AP/ vs. log (nj no), including all data from panel (b). 

— firiQ 1 exp(AP//3) or, with a cutoff, dAP/dn = 
-/3n " 1 [exp(AP//3) - 1], which prod uces the logarithmic 
form of Eq. [2] for small n, with saturation at n = uq. 
This suggests an activated process, perhaps involving 
a generalized ensemble, such as the stress ensemble, as 
discussed by several authors p~7H2Q] . 

To summarize: First, for frictional granular systems 
in/near the shear jamming regime, </>s < </> < </>j, we 
generated sheared states without shear bands, even with 
large strains or over many cycles of shear, making it 
possible to experimentally probe the constitutive rela- 
tions of granular materials. Second, we have character- 
ized a Reynolds effect for fixed (j) that is approximately 
quadratic in 7 using R = (d 2 P / d^) /2, and R has a 
power-law scaling with A<j). We note that the specific 
form for P(j) may well depend on the particle interac- 
tion force; a more general form might be P = R(A(j) a )j 5 , 
where for our experiments, a ~ —3.3 and 5 ~ 2. Third, 
we found that under cyclic shear, the stresses in fric- 
tional granular systems evolve logarithmically slowly, as 
one might expect for an activated process, towards a 
state where the pressure is symmetric, modulo fluctua- 
tions, about the mid-point of strain. The pressure at the 
symmetry point may not be zero. Fourth, this slow evo- 
lution is characterized by highly novel scaling behavior, 
such that there is good collapse of all data. 

These results point towards several interesting direc- 
tions. First, it is reasonable to search for a description of 
these states in terms of an ensemble picture, such as the 
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